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1. FORMAT 

Text should be produced within the dimensions 
shown on these pages: each column 7.5 cm wide 
with 1 cm middle margin, total width of 16 cm 
and a maximum length of 20.2 cm on first pages 
and 21 cm on second and following pages. The 
ETpX document style uses the maximal stipu- 
lated length apart from the following two excep- 
tions (i) DTpX does not begin a new section di- 
rectly at the bottom of a page, but transfers the 
heading to the top of the next page; (ii) MpX 
never (well, hardly ever) exceeds the length of the 
text area in order to complete a section of text or 
a paragraph. 

1.1. Spacing 

We normally recommend the use of 1.0 (sin- 
gle) line spacing. However, when typing com- 
plicated mathematical text MpX automatically 
increases the space between text lines in order 
to prevent sub- and superscript fonts overlapping 
one another and making your printed matter il- 
legible. 

1.2. Fonts 

These instructions have been produced using 
a 10 point Computer Modern Roman. Other 
recommended fonts are 10 point Times Roman, 



'Footnotes should appear on the first page only to indi- 
cate your present address (if different from your normal 
address), research grant, sponsoring agency, etc. These 
are obtained with the \thanks command. 



New Century Schoolbook, Bookman Light and 
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2. PRINTOUT 

The most suitable printer is a laser printer. A 
dot matrix printer should only be used if it pos- 
sesses an 18 or 24 pin printhead ("letter-quality"). 

The printout submitted should be an original; 
a photocopy is not acceptable. Please make use 
of good quality plain white A4 (or US Letter) 
paper size. The dimensions shown here should be 
strictly adhered to: do not make changes to these 
dimensions, which are determined by the docu- 
ment style. The document style leaves at least 
3 cm at the top of the page before the head, which 
contains the page number. 

Printers sometimes produce text which con- 
tains light and dark streaks, or has considerable 
lighting variation either between left-hand and 
right-hand margins or between text heads and 
bottoms. To achieve optimal reproduction qual- 
ity, the contrast of text lettering must be uniform, 
sharp and dark over the whole page and through- 
out the article. 

If corrections are made to the text, print com- 
pletely new replacement pages. The contrast on 
these pages should be consistent with the rest 
of the paper as should text dimensions and font 
sizes. 
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Table 1 



Biologically treated effluents (mg/1) 




Pilot plant 




Full scale plant 






Influent 


Effluent 


Influent 


Effluent 


Total cyanide 


6.5 


0.35 


2.0 


0.30 


Method-C cyanide 


4.1 


0.05 




0.02 


Thiocyanide 


60.0 


1.0 


50.0 


< 0.10 


Ammonia 


6.0 


0.50 




0.10 


Copper 


1.0 


0.04 


1.0 


0.05 


Suspended solids 








< 10.0 


Reprinted from: G.M. 


Ritcey, Tailings 


Management, Elsevier, 


Amsterdam, 





1989, p. 635. 



3. TABLES AND ILLUSTRATIONS 

Tables should be made with I^TjrjX; illustra- 
tions should be originals or sharp prints. They 
should be arranged throughout the text and 
preferably be included on the same page as they 
are first discussed. They should have a self- 
contained caption and be positioned in flush- 
left alignment with the text margin within the 
column. If they do not fit into one column 
they may be placed across both columns (us- 
ing \begin{table*} or \begin{f igure*} so that 
they appear at the top of a page). 

3.1. Tables 

Tables should be presented in the form shown 
in Table |l|. Their layout should be consistent 
throughout. 

Horizontal lines should be placed above and be- 
low table headings, above the subheadings and at 
the end of the table above any notes. Vertical 
lines should be avoided. 

If a table is too long to fit onto one page, 
the table number and headings should be re- 
peated above the continuation of the table. For 
this you have to reset the table counter with 
\addtocounter-[table}{-l}. Alternatively, the 
table can be turned by 90° ('landscape mode') 
and spread over two consecutive pages (first 
an even-numbered, then an odd-numbered one) 
created by means of \begin{table> [h] without 
a caption. To do this, you prepare the table as 
a separate ETgX document and attach the tables 
to the empty pages with a few spots of suitable 
glue. 



3.2. Line drawings 

Line drawings should be drawn in India ink 
on tracing paper with the aid of a stencil or 
should be glossy prints of the same; computer pre- 
pared drawings are also acceptable. They should 
be attached to your manuscript page, correctly 
aligned, using suitable glue and not transparent 
tape. When placing a figure at the top of a page, 
the top of the figure should be at the same level 
as the bottom of the first text line. 

All notations and lettering should be no less 
than 2 mm high. The use of heavy black, bold 
lettering should be avoided as this will look un- 
pleasantly dark when printed. 

3.3. Black and white photographs 

Photographs must always be sharp originals 
(not screened versions) and rich in contrast. 
They will undergo the same reduction as the text 
and should be pasted on your page in the same 
way as line drawings. 

3.4. Colour photographs 

Sharp originals (not transparencies or slides) 
should be submitted close to the size expected 
in publication. Charges for the processing and 
printing of colour will be passed on to the 
author(s) of the paper. As costs involved are per 
page, care should be taken in the selection of size 
and shape so that two or more illustrations may 
be fitted together on one page. Please contact 
the Technical Editor in the Camera-Ready Pub- 
lications Department at Elsevier for a price quo- 
tation and layout instructions before producing 
your paper in its final form. 
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Figure 1. Good sharp prints should be used and 
not (distorted) photocopies. 



in parentheses flush with the right-hand margin 
of the text and level with the last line of the equa- 
tion. For multi-line equations, use the eqnarray 
environment. For complex mathematics, use the 
A\4S-WT^i package. 
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References should be collected at the end of 
your paper. Do not begin them on a new page 
unless this is absolutely necessary. They should 
be prepared according to the sequential numeric 
system making sure that all material mentioned 
is generally available to the reader. Use \cite 
to refer to the entries in the bibliography so that 
your accumulated list corresponds to the citations 
made in the text body. 

Above we have listed some references according 
to the sequential numeric system Jl]-||] . 



Figure 2. Remember to keep details clear and 
large enough. 



4. EQUATIONS 

Equations should be flush-left with the text 
margin; DTgX ensures that the equation is pre- 
ceded and followed by one line of white space. 
DT^X provides the document-style option f leqn 
to get the flush- left effect. 

H aP {u) = E^(co)S aP + (a\W % \p) (1) 

You need not put in equation numbers, since 
this is taken care of automatically. The equation 
numbers are always consecutive and are printed 
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Vector and Axial- Vector Spectral Functions and QCD 

Andreas Hocker a> 

aa Laboratoire de l'Accelerateur Lineaire, 

IN2P3-CNRS et Universite de Paris-Sud, F-91405 Orsay, France 

We present new results for the r hadronic spectral functions analysis using data accumulated by the ALEPH 
detector at LEP during the years 1991-94. In addition to the vector spectral function, the axial-vector spectral 
function and, separately, the r~ — » 2tt~tv + v t as well as the r~ —* iv~2n v T spectral functions are determined from 
their respective unfolded, i.e., physical invariant mass spectra. The spectral functions are applied to QCD chiral 
sum rules in order to extract information about saturation at the r mass scale. Using the the semi-leptonic r 
decay rate for vector and axial-vector currents in addition to spectral moments, we obtain precise measurements 
of the strong coupling constant a s (M T ) and the contributing non-perturbative power terms. The evolution to the 
Z mass yields a s (M z ) = 0.1219 ± 0.0019. 



1. INTRODUCTION 

Spectral functions of hadronic r decays are the 
key objects for various interesting studies con- 
cerning resonance structure analysis, the Con- 
served Vector Current hypothesis (CVC) and 
QCD tests involving strong sum rules and the 
measurement of a s . The ALEPH measurement 
of the vector spectral functions and related topics 
was presented at the same conference by R. Ale- 
many [Q (see also 

In this article we deal with new data of the 
spectral functions from non-strange, axial-vector 
hadronic r decays measured by the ALEPH Col- 
laboration. We create the sum and the difference 
of vector and axial-vector spectral functions to 
get access to information about the saturation of 
QCD sum rules at the r mass scale. As an ap- 
plication, we determine the electric polarisability 
of the pion from chiral QCD constraints. In cal- 
culating the semi-leptonic r widths and so-called 
spectral moments of the vector and axial-vector 
currents, we independently fit the strong cou- 
pling constant a s (M T ) and the respective non- 
perturbative contributions of the Operator Prod- 
uct Expansion (OPE) power series |3|,|| using the 
framework already exploited in previous analy- 
ses |3]-fj]]. We finally perform a combined fit of 
all components in order to obtain the best value 
of a s {M T ). 



2. SPECTRAL FUNCTIONS OF VEC- 
TOR AND AXIAL-VECTOR CUR- 
RENTS 

The measurement of the non-strange r vec- 
tor (axial-vector) current spectral functions re- 
quires the selection and identification of r de- 
cay modes with a G-parity G=+l (G=— 1), i.e., 
hadronic channels with an even (odd) number of 
neutral or charged pions. The isovector spectral 
function Vi y- (a% a-) of a vector (axial- vector) 
r decay channel V~ v T (A~ v T ) is obtained by 
dividing the normalized invariant mass-squared 
distribution (l/N v / A - )(dN v / A - /ds) for a given 
hadronic mass \fs by the appropriate kinematic 
factor (throughout this article, charge conjugate 
states are implied): 



«l,V-Ah,yl-(s) = 

M 2 T B(t- 



V/A~ v r ) 



6|Kd| 2 S EW B(t- 



V/A- 



N- 



V/A' 



ds 



e v T v e ) 

— Y (\ — 



-1 
,(1) 



where \V ud \ = 0.9752 ± 0.0007 denotes the 
CKM weak mixing matrix element j|] an d 
Sew — 1-0194 accounts for electroweak sec- 
ond order corrections The r mass M T = 



1776.961^7 MeV/c 2 is 
cent BES measurement 



taken 



from the re- 
The spectral 
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functions are normalized by the ratio of the 
respective vector/axial- vector branching frac- 
tion B(t~ — > V~ l 'A" v T ) to the branch- 
ing fraction of the electron channel B(r~ — ► 
e~D e v T ) = 17.79 ±0.04 Q, where the latter is 
additionally constrained via universality from 
B(t~ — > and the r lifetime. Note that 

our definition of the r spectral functions differs 
from the one used in jl2j by an additional factor 
4tt 2 . 

Assuming unitarity (which implies the opti- 
cal theorem) and analyticity to hold, the spec- 
tral functions of hadronic r decays are related 
via dispersion relations to the imaginary parts 
of the two-point correlation functions 11^- v (q) = 
i J d^x e lqx (0\T(U^(x)U»(0)^)\0) of vector (U = 
V = V-'j7 A, '0i) or axial-vector (U = A = 
■0^7^75 V'i) colour-singlet quark currents in cor- 
responding quantum states (see, e.g., PJl^|). 

3. THE MEASUREMENT PROCEDURE 

The measurement of the spectral functions 
defined in Eq. (|l]) requires the determination of 
the physical invariant mass-squared distribution. 
The details of the analysis are reported in 0. A 
description of the ALEPH detector and its per- 
formance is published in [Q. 

In the following, we will briefly outline the im- 
portant steps of the measurement procedure: 

- Tau pairs originating from Z° decays are 
detected utilizing their characteristic collinear jet 
signature and the low multiplicity of their de- 
cays. Using the data from 1991-94, a total 
124 358 r pairs is selected corresponding to a de- 
tection efficiency of (78.8 ±0.1)%. The overall 
non-r background contribution in the hadronic 
modes amounts to (0.6 ±0.2)%. Details about 
the ALEPH r pair selection are given in [l^-]T8| . 

- Charged particles (electrons, muons and 
hadrons) are identified employing a maximum 
likelihood method to combine different and es- 
sentially uncorrelated information measured for 
each individual track. The procedure and the 
discriminating variables used in this analysis are 
described in Jl^Jl^]. The calibration of charged 



tracks is performed using low radiating e + e~— ► 
H + ji~ events at beam energy and the invariant 
mass measurement of well-known, narrow reso- 
nances at low and intermediate energies. The 
resulting calibration uncertainty amounts to less 
than 0.1%. 

Photons are reconstructed by collecting 
associated energetic electromagnetic calorimeter 
(ECAL) towers, forming a cluster. To distinguish 
genuine photons from fake photons a likelihood 
method is applied using ECAL information, e.g., 
the fraction of energy in the respective ECAL 
stacks, the transverse size of the shower or the 
distance between the barycentre of the cluster 
and the closest charged track. The energy cali- 
bration is performed using electrons originating 
from Bhabha, r and two-photon events. We ob- 
tain a relative calibration uncertainty of about 
1.5% at low energy, 1% at intermediate energies 
and 0.5% at high energy. 

- The 7T° finder uses the virtue of a 7r°-mass 
constraint fit to correctly attribute two recon- 
structed photons to the corresponding 7r° decay. 
At higher 7r° energy, the opening angle between 
the boosted photons tends to become smaller 
than the calorimeter resolution so that the two 
electromagnetic showers are often merged in one 
cluster. The transverse energy distribution in 
the ECAL nevertheless allows the computation of 
energy- weighted moments providing a measure of 
the two-photon invariant mass. Remaining pho- 
tons are considered as originating from a ir° where 
the second photon has been lost. 

- The classification of the inclusive hadronic 
r decay channels is performed according to on 
the basis of the number of reconstructed charged 
and neutral pions. The exclusive channels listed 
in Table [l] are obtained by subtracting the r and 
non-r background and the strange contribution 
from the inclusive measurements using the Monte 
Carlo simulation. In order to extract the physical 
invariant mass spectra from the measured ones it 
needs to be unfolded from the effects of measure- 
ment distortion. 

- The unfolding method used here is based 
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Table 1 

Vector and axial-vector hadronic r decay modes with their contributing branching 
fractions. According to jL5|, the KKtt channels are assumed to be (78l2s)% vec tor 
and (25 ± 25)% axial-vector, while the errors are 100% anti-correlated. The KKtttt v t 
modes are conservatively assumed to have (50 ± 50)% vector and axial-vector parts. 



Vector 


BR (in %) 


Axial- Vector 


BR (in 7o ) 


_ n 
7T 7T V T 


25.35 ±0.19 




11.23±0.16 


7T~37r° V T 


1.17±0.14 


tt-2tt v t 


9.23±0.17 


27r~7r + 7r° v t 


2.54 ±0.09 


2tt~it + v t 


9.13 ±0.15 


TT~hlT {) V T 




7r~47T° V T 


0.03±0.03( 1 ) 


2tt~tt+3tt°v t 


| 0.04 ± 0.02^) 


2-K-TT+21T V T 


0.10±0.02 


37T-27r+7r v T 




37T~27T + V r 


0.08 ±0.02 


UJ1T~ V T ^ 


1.83 ±0.09 


ujir~ir° v T ^ 


0.42 ±0.09 


?7 7r~7r v T ^ 


0.17±0.03 


Tj 27T~7T + V T 


0.04 ±0.01 






rj 7r"27r° v T 


0.02 ±0.01 


K- K u v T 


0.19 ±0.04 








0.13±0.05 


KTK+TT- V T 


0.04 ±0.04 


K°K 7T- V T 


0.13±0.05( 1 ) 


K°K 7T- V T 


0.04 ±0.04 


KK°7T V T 


0.08 ±0.04 


KK°7r° v T 


0.03 ±0.03 


KK7T7T 


0.08±0.08 (1) 


KK7T7T 


0.08±0.08 (1) 


Total Vector 


31.71 ±0.31 


Total Axial- Vector 


30.42 ±0.32 



1 The branching ratio is obtained using constraints from isospin symmetry (see text and |2| ) . 

2 Through w — » 7r~7r+7r , 88.8% of this channel is reconstructed in 2tt~ 7r+7r° u T . 

3 Through r\ — •> 2"j, 39.3% of this channel is reconstructed in tt~ 37r° v T - 



on the regularised inversion of the detector re- 
sponse matrix, obtained from the Monte Carlo 
simulation, using the Singular Value Composition 
technique. Details are published in p(| . 

Systematic Errors 

The study of systematic errors affecting the 
measurement is subdivided into several classes ac- 
cording to their origin, viz., the photon and 7T° 
reconstruction, the charged track measurement, 
the unfolding procedure and additional sources. 
Since we use an unfolding procedure based upon 
a detector response matrix from the Monte Carlo 
simulation, the reliability of the simulation has to 
be subjected to detailed studies p[,p]JTq]. 

In order to check the photon reconstruction in 
the ECAL, we studied the influence of calibra- 
tion and resolution uncertainties as well as pos- 
sible problems with the reference distributions of 



a likelihood used to veto fake photon candidates, 
of the contamination and energy distribution of 
fake photons, of the photon detection efficiency 
at threshold energies (£t hresh = 300 MeV) and 
in the neighborhood of charged tracks. 

Correspondingly, the effects of calibration and 
resolution uncertainties in the measurement of 
charged tracks are checked, accompanied by tests 
of the reconstruction efficiency of highly colli- 
mated multi-prong events, and the simulation of 
secondary nuclear interactions. 

In addition, systematic errors introduced by 
the unfolding procedure are tested by compar- 
ing known, true distributions to their respective 
unfolded ones and by varying the regularisation 
conditions. 

Finally we add systematic errors due to the 
limited Monte Carlo statistics and to uncertain- 
ties in the branching ratios of the respective r 
decay modes. 
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Invariant Mass Spectra and Spectral Func- 
tions 

The measurement procedure sketched above 
provides the physical invariant mass spectra of 
the measured r decay modes including their re- 
spective bin-to-bin covariance matrices obtained 
after the unfolding of the spectra from the sta- 
tistical errors and the study of systematic uncer- 
tainties. 

The exclusive vector and axial-vector r de- 
cay channels are listed in Table [l| The branch- 
ing ratios shown are mainly taken from Jl8| , in- 
cluding recent improvements in the precision of 
the measurements of r decays involving r\ and u> 
mesons 21 2^] as well as Kaons p3|. 

The total vector and axial-vector current spec- 
tral function (|]J) is obtained by summing up 
the exclusive spectral functions with the ad- 
dition of small contributions from unmeasured 
modes. These are taken from the K0RALZ3 . 8 , 
TAU0LA1 . 5 |24| r Monte Carlo simulation (ac- 
companied by an inflated systematic error), as 
discussed below. 

Vector Spectral Function 

The two- and four-pion states are measured ex- 
clusively, while the six-pion state is only partly 
measured. The total six-pion branching ratio has 
been determined in Q using isospin symmetry. 
However, this procedure is not straight forward, 
as one has to consider that the six-pion channel is 
contaminated by isospin- violating t~ — > r\ 2tt~tt + , 
7/7r~27r°i/ T decays (which were reported for the 
first time by the CLEO Collaboration (22)). 

The small fraction of the to 7r~ v T decay channel 
that is not reconstructed in the four-pion final 
state is added using the simulation. 

Similarly, we correct for r\ ir~ir° v T decay modes 
other than r\ — > 27 which is classified in the 



h-37T U 



final state. 



The K~ K° v T is taken entirely from the simu- 
lation. 

Using isospin constraints from r branching ra- 
tios deduced in (25), it can be shown yl| that 
the KK-7T channels are (78jl||)% vector. The cor- 
responding spectral functions are obtained from 
the Monte Carlo simulation. 



. if ^(VM=l)v T (ALEPH) 
W eV -4 (V, 1=1) 

naive parton model prediction 




Mass 2 (GeV/c 2 ) 2 



Figure 1. Total r vector current spectral func- 
tion. The band shows the corresponding e + e~ 
annihilation data obtained from isovector final 
states translated to r decays via isospin rotation 
as comprehensively explained in The dashed 
line represents the naive parton model prediction, 
while the QCD corrected prediction lies roughly 
20% higher. 



A preliminary ALEPH analysis of r decays 
into KK7T7T final states yields a total branching ra- 
tio of (0.17 ± 0.05)%. However, their vector and 
axial- vector parts are unknown and are therefore 
both estimated to be (50 ±50)%. 

For additional details on the vector spectral 
functions see Ref. The total r vector 

spectral function is compared to the isospin ro- 
tated cross sections from isovector e + e~ final 
states ||[[| is depicted in Fig. |l|. The agreement 
between both curves is fairly good, keeping in 
mind the appearance of large bin-to-bin correla- 
tions within the r spectral function as a conse- 
quence of the necessary regularisation during the 
unfolding procedure. 
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Figure 2. Unfolded (physical) invariant mass- 
squared spectra of the r final states 2tt^tt + v t 
and 7T~27r° f r . The dashed/dotted line shows the 
corresponding measured distribution. 



Figure 3. Total r axial- vector current spec- 
tral function. The dashed line represents the 
naive parton model prediction, while the QCD- 
corrected prediction lies roughly 20% higher. 



Axial- Vector Spectral Function 

The three-pion spectral functions are measured 
in both occurring final states. Fig. ^ shows 
the unfolded 27r~7r + v t and 7r~27r° v T mass spec- 
tra with reasonable agreement in form and nor- 
malization. The small bump of the measured 
27r~7r + v T spectrum at low mass (Fig. |J) is caused 
by decays where only two charged tracks are re- 
constructed. Due to incomplete ECAL energy 
collection, the measured 7r - 27r°zy T distribution is 
slightly shifted to lower masses. These details are 
well reproduced by the Monte Carlo simulation. 
We assume in the following that both channels 
have identical spectra and use the weighted aver- 
age of the distributions for the total axial-vector 
spectral function. 

The five-pion spectral functions are exclusively 
measured in the 2tt~tt + 2tt° iy T and 3ir~2ir + v T 
final states. Using Pais' isospin classes ^6|, the 
branching fraction of 7r~47r° v T can be bounded 
entirely using the 37r _ 27r + 7r° v T branching frac- 
tion. It is found to be smaller than 0.055%. We 



take half of this upper limit with an error of 100%. 

In analogy to the vector case, we add the small 
fraction of the ujir~ir v T decay channel that is 
not accumulated in the 27r _ 7r + 27r° v T final state 
using the simulation. 

Also considered are the axial-vector ?7(37r)~ v T 
final states |%| and, as in the vector case, a 
(50 ±50)% KKtttt contribution. Both spectral 
functions are taken from the simulation accom- 
panied by comfortable systematic errors due to 
the uncertainty of the respective invariant mass 
distributions. 

The total axial- vector spectral function is plot- 
ted in Fig. [|. However, we are still not in the 
asymptotic region at the r mass scale. We may 
expect additional oscillations to lift the spectral 
function to roughly 1.2 times the naive parton 
model prediction (that is expected asymptoti- 
cally). 



G 




Figure 4. Inclusively measured vector plus axial- 
vector (vi + <zi) spectral function. Again, the 
dashed line represents the naive parton model 
prediction. 



(i>i + ax) Spectral Function 

In the favourable case of the vector plus axial- 
vector spectral function we do not have to dis- 
tinguish the current properties of the respective 
non-strange hadronic t decay channels. Hence 
we measure the mixture of all contributing non- 
strange final states as inclusively as possible. The 
dominating two- and three-pion final states are 
still measured exclusively, while the remaining 
contributing topologies are treated inclusively, 
i.e., without any subtraction of r-background 
originating from one of the considered decay 
modes. This improves the statistical error. An- 
other advantage is, that we do not have to worry 
about the current properties of the KK-7T and 
KK7T7T modes or about possible missing, i.e., un- 
measured, t decay modes as they are necessarily 
contained in the inclusive sample. 

The + ai) spectral function is depicted in 
Fig. ^. The improvement in precision in com- 
parison to an exclusive sum of Fig. and Fig. | 



Figure 5. Vector minus axial- vector (vi — a{) 
spectral function. The band shows the corre- 
sponding spectral function using two times the 
vector spectral function from e + e~ annihilation 
data minus the inclusively measured (i>i + a{) 
spectral function from r decays. 



becomes obvious at higher mass-squared. One 
clearly sees the oscillating spectral function which 
does not seem to approach the asymptotic limit 
at s — > Af^, which is predicted from perturba- 
tive QCD to lie about 20% higher than the naive 
parton model prediction. 



4. QCD CHIRAL SUM RULES 

As important application of the results ob- 
tained from the spectral function analysis above, 
we present a study of four QCD sum rules in- 
volving the difference of the vector and the axial- 
vector spectral functions, each deduced from vec- 
tor and axial-vector current conservation in the 
limit of the chiral symmetry (m u — = 0, 
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m 



47T 



' Z' |fc i;i(a)-oi(a) = ^ 2 



o 

so — >°o 



4tt 2 



ds(vi{s) - ai(s)) = 



- F A , (2) 



(3) 



Sq— »OC 

^2 y dss (t>i(s) - ai(s)) = 0, 
o 

so — >oo 

— / dssln-p (vi(s) -ai(s)) = 
16tt 2 / 2 



(4) 



3 a 



(5) 



Eq. (g) is known as the Das-Mathur-Okubo 
(DMO) sum rule |27|. It relates the given inte- 
gral (Idmo) to the square of the pion decay con- 
stant!] / 7r = (92.4 ±0.26) MeV @ obtained from 



the decays n 



and 7r 



H VyTi\ to 



the pion axial-vector form factor Fa for radia- 
tive decays tt~ — > l~vci\ and to the pion charge 
radius-squared (r 2 ) = (0.431 ± 0.026) fm 2 {29), 
determined by means of a careful fit to space- 
like data [^8) using a two-loop chiral expan- 
sion coefficient as additional degree of free- 
dom. The error of (r 2 ) includes theoretical un- 
certainties. For completeness we will also ex- 
press the results on the pion polarisability (see 
next paragraph) in terms of the standard value 
(rl) = (0.439 ±0.008) fm 2 @ obtained from a 
one parameter fit to the same data ^8|. Eqs (||) 
and (|J) are the first and second Weinberg sum 
rules (WSR) j30), where the pion pole in (||) is 
already integrated out. When switching quark 
masses on, only the first WSR remains valid while 
the second WSR breaks down due to contribu- 
tions from the difference of non-conserved vec- 
tor and axial- vector currents of order m^/s, lead- 
ing to a quadratic divergence of the integral . 

1 We now identify the vector and axial-vector spectral 
functions with the absorptive parts of the corresponding 
two-point correlation functions II'"'.,, . . 

* ' ri 
Our definition of /V differs from the one used in [pj by 

a factor of \/2- 



Finally, Eq. (|J) represents the electromagnetic 
splitting of the pion masses |3^| , 

Fig. U shows the vector minus axial- vector dis- 
tribution obtained from t decays in addition to 
two times the vector spectral function from e + e~ 
annihilation data minus the inclusive (vi + a±) 
measurement, yielding a corresponding (vi ~ ai) 
spectral function with a slightly better preci- 
sion at the end of the r phase space. The sum 
rules (§)-(|^) versus the upper integration bound 
so < Mj? are plotted in Figs. ^a-d. In addition to 
the the r spectral functions, we again show the 
contribution from e + e~ annihilation combined 
with the (v\ + a\) inclusive measurement. The 
horizontal band depicts the corresponding chiral 
prediction of the integrals taken from []l2| . 

We stress that only for the DMO sum rule 
(Fig. ^a), where contributions from higher mass- 
squares are suppressed, does the saturation 
within the one sigma error seem to occur at the 
r mass scale. The other sum rules (Fig. ||b-c) 
are apparently not saturated at M 2 as already 
indicated by the non- vanishing (vi — ai) spectral 
function at the end of the r phase space. This 
is due to a small axial- vector contribution that is 
clearly below the QCD prediction at M 2 as can 
be seen in Fig. S. 

The Polarisability of the Pion 

In order to quantify the actual precision of the 
sum rules above, we determine the electric polar- 
isability^ of the pion, given by J|3| 



ocFa 



(0) 



utilizing the DMO sum rule Idmo @) a $ 
proposed in 34 1. The computation of the 
pion axial-vector form factor Fa is most ac- 
curate using the measurement of the ratio 
7 = Fa/Fv obtained from radiative pion de- 
cays 7r~ —>■ e~v e j. Using the weighted aver- 
age 7 = 0.46 ±0.05 of the measurements 35 



3 The electric polarisability c?e of a physical system can 
be understood classically as the proportionality constant 
which governs the induction of the dipole moment p of 
a system in the presence of an external electric field E: 
p = «eE. The polarisability is an important quantity to 
characterize a particle, i.e., in probing its inner structure. 
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Figure 6. Sum rules corresponding to Eqs (|^)-(^) (plots: a-d) versus the upper integration bound sq. 



as well as the CVC relation between the pion 
vector form factor and the 7r° lifetime \Fy\ = 
(1/a) (2itt„o m n o )" 1 / 2 = 0.0132 ± 0.0005 Ql , 
we have 

af eo = (2.86 ± 0.33) x 10~ 4 fm 3 (7) 

as a theoretical prediction. 

Using the r (v% — ai) spectral functions, we 
evaluate the DMO integral at the r mass scale as 
I£ uo = (28.0 ± 1.6 ± 1.1) xlO" 3 , where the first 
error is the direct propagated integration er- 
ror using the covariance matrices of the spec- 
tral functions. The second error accounts for 
interpolation biases that occur because we do 
not integrate up to exactly sq = M 2 since 



the error of the spectral function diverges at 
the very end of the r phase space. For the 
2«i(e + e~)— (v\ + ai)(r) spectral function we 

obtain 7d'm + o _ =(25.2 ±1.3 ±0.4) xlO" 3 with a 
X 2 = 2.7/1 when compared to the above pure 
t result, assuming both measurements to be 50% 
correlated. The weighted average (taking into ac- 
count the correlation) is 

(/dmo) = (26.4 ± 1.5) x 10- 3 . (8) 

According to the prescription of Ref. ||, we in- 
creased the error by X 2 /l to account for some 
inconsistency. Combining (||) and (ph with the 
assumption that the contribution to (H) for sq > 
M 2 is negligible, i.e., the integral is saturated, 
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one finds that the pion polarisability is 

= (2.40 [2.68] ± 1.14 [0.76]) x 10~ 4 fm 3 . (9) 

The figures in brackets give the correspond- 
ing result if we use the standard value of 
(r 2 ) = (0.439 ± 0.008) fm 2 |§ for the pion charge 
radius-squared. Both results @ are in agreement 
with the chiral prediction ((?]). 

The authors of |m[ (see also |37j) used the first 
WSR © as an additional constraint to consider- 
ably improve the precision and the reliability of 
the /dmo evaluation as it naturally reduces the 
sensitivity to the saturation assumption. 

Another approach to the solution of (|J) deals 
with the Laplace-transformed DMO sum rule, 
which suppresses the high energy tail of the spec- 
tral function in order to improve saturation at 
M 2 and to increase the precision of the inte- 
gral ||||: 



s/M 2 Ms) - ai(s) 
s 



fl C 6 (0(6)) C 8 (0(8)) 



M 2 6M 6 



24M 8 



(10) 



where /dmo 

(M 2 ) = /dmo in the limit M — ► cx>. 
At sufficiently high M 2 the impact of the di- 
mension D = 6 and dimension D = 8 non- 
perturbative terms on 7dmo(-^ 2 ) is small. One 
may use projecting sum rules with improved 
saturation to determine the corresponding phe- 
nomenological operators: 

so 

-47r 2 C 6 (0(6)) ~ Jdss 2 {v 1 {s)-a 1 (s)) 
o 

so 

- Pisojdss (vx(s) - ai(s)) , (11) 





47r 2 C 8 (C(8)) ~ / ds s">i(s) - 01(a)) 
o 

so 

-P 2 s^Jdss(v 1 (s)-a 1 (s)) . (12) 



Obviously, Eq. (|ll|) and (gj) hold for s — > oo 
since the second term on the r.h.s. vanishes in 
the chiral limit by virtue of the second WSR (||) . 
Unfortunately, the gain in convergence of the sum 
rules obtained from the insertion of the (even less 
convergent) second WSR is not very successful 
as it is accompanied by large additional errors. 
The coefficients (3i/ 2 in and ( p"2| ) depend 

on the high energy tail of the (v\ — ax) spec- 
tral function which can be expanded in powers 
of so- They are estimated to be f3\ ~ 1.5 and 
/?2 ~ 2.8. We obtain, for the non-perturbative 
contributions, the extremely rough estimates (us- 
ing t data only): C 6 (O(6))=0.025 ± 0.027 GeV 6 
and C 8 (0(8))=-O.15±O.16 GeV 8 . This yields 
the most precise value of Idmo=(26.3 ± 0.6) for 
r spectral functions at (high) M 2 =2.6 GeV 2 en- 
abling one to avoid the large errors of the inaccu- 
rate operator terms at low energy. 



A more promising solution of Eq. (10) lies 
in a simultaneous fit of /dmo and the non- 
perturbative terms by means of moments on the 
Borel parameter A/ 2 . We perform a x 2 fit using 
the (strongly correlated) moments M 2 = 0.2, 0.7, 
1.2, . . ., 3.7 GeV 2 in order to guarantee sufficient 
information for the constraint of the dimension 
D = 6 and D — 8 operators at low M 2 and non- 
biased, purely perturbative contributions at high 
M 2 to fix /dmo through the Laplace-transformed 
DMO integral in @. The fit converges with a 
X 2 = 2.3 over 5 degrees of freedom, yielding (for 
r data) 



/dmo = (25.8 ± 0.3 ± 0.1) x 10 



-3 



(13) 



which is in agreement with (g). The sec- 
ond error accounts for estimated uncertain- 
ties induced by the saturation assumption. 
From the fit, the non-perturbative contribu- 
tions are C 6 (C(6)) =0.0029 ±0.0002 GeV 6 and 
C 8 (O(8))=-0.0015± 0.0003 GeV 8 with an anti- 
correlation of nearly 100%. 

Fig. depicts /dmo as a function of M 2 for the 
pure r and the combined r and e + e~ annihilation 
(vi — a%) spectral functions. In this case, the pure 
r spectral function provides smaller errors be- 
cause of the more accurate two-pion contribution 
in r decays which is strongly weighted in (fl0|). 
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Figure 7. Laplace-transformed DMO sum rule 
Idmo as a function of the Borel parameter M 2 for 
the exclusive r and the combined e + e~, r <zi) 
spectral function. The shaded area gives the r 
integral without the non-perturbative contribu- 
tions. It diverges hyperbolically at small M 2 . 



The shaded band shows the r result without the 
non-perturbative terms. Non-perturbative con- 
tributions dominate the integral at M 2 < 1 GeV 2 
and become negligible for M 2 > 2.5 GeV 2 as the 
comparison between the shaded and the hatched 
bands demonstrates. 

Using ( |l3| ) the electric polarisability is 

a c E xp = (2.68 [2.96] ± 0.91 [0.32]) x 10~ 4 fm 3 (14) 

with significantly lower errors compared to @ 
due to the improved convergence of the Laplace- 
transformed DMO sum rule. Again the figures in 
brackets give the corresponding result when using 
the standard value of (r 2 ). 

5. THE MEASUREMENT OF a s 

Tests of Quantum Chromodynamics and the 
precise measurement of the strong coupling con- 
stant a s at the r mass scale carried out for the 
first time by ALEPH §§ and CLEO § have 



been the subject of engaged discussions about 
theoretical and experimental implications, ac- 
companied by a considerable number of inter- 
esting publications (see, e.g., |3l}fl0| . |l^ , |4lf ) . Fol- 
lowing the spirit of the previous ALEPH analy- 
sis [gjj^], we determine a s (M T ) and the contribut- 
ing non-perturbative terms for the first time us- 
ing the vector, axial- vector as well as vector plus 
axial-vector spectral functions and from a com- 
bined fit of all substantially uncorrelated input 
variables. 

The t Semileptonic Width 

Using the analytic property of the vector/axial- 
vector two-point correlation functions Tl^ v ^ A , 
the non-strange, semileptonic r decay rates 
R T y/A can be expressed as contour integrals in 
the complex s-plane [fl2"|.[H 



Rt,V/A 



6ni 



ds 
M 2 



2s \ tt(0+1) 
M 2 J vd,V/A 



-Y 

M 2 ) 



2s tt(0) 

j^j2 ud,A 



(15) 



The advantage of this formulation is that it avoids 
the integration over the whole energy-squared 
spectrum which includes large, uncontrolled non- 
perturbative effects. 

By means of the short distance OPE, the the- 
oretical prediction of R t _v/a separates in the fol- 
lowing contributions 

3 

Rt,v/a = 2 \Vud\ Sew(1 + <5ew + 

Sp + tfmass + <5np) , (16) 

where \V ud \ = 0.9752 ± 0.0007 [|, Sew = 
1.0194 as in pi) and the small non-logarithmic 



electroweak correction 8' % 



EW 



0.0010 



The perturbative contribution <5p of the corre- 
lators in (|l5| ) is known to next-to-next-to-leading 



order 



Moreover, the series is resummcd 



to all orders up to the unknown perturbative 
coefficients of the QCD /3-function. Hence the 
convergence of the series is reinforced and the de- 
pendence on the scale and the choice of the renor- 
malization scheme is reduced [Eof . 
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The leading quark mass correction 5 mass in (16) 
is safely neglected for u and d quarks (« 
-0.08% @). 

For the non-perturbative contribution <5np 
we use the SVZ-approach H in terms of 
a power series of M T . Long-distance QCD 
effects are thereby separated into unpredictable, 
non-perturbative contributions, absorbed in di- 
mensional, phenomenological operators O(D) 
and into short distance parts Cd (Wilson 
Coefficients), calculable within perturbative 
QCD: 



= E siD) = E n 



0(D) 



D=4,6,8 



Z>=4,6,8 



(„ M 2)D/2 



(17) 



The first operator that involves non-perturbative 
physics appears at D = 4 and is linked to the 
gluon condensate. 

Spectral Moments 

To detach the measurement from theoretical 
constraints on the operators O(D), it is conve- 
nient to fit simultaneously a s (M T ), 0(4), 0(6) 
and 0(8). In order to incorporate new experi- 
mental information, we exploit the explicit shape 
of the (normalized) vector/axial-vector invariant 
mass spectrum (1/iVy/^) dNy/^/ds by means of 
the spectral moments 

Ml k l 

(£) 



1 dWi 



V/A 



VI A 



ds 



(18) 



where the choice of k = 1 and I = 0, 1, 2, 3 pro- 
vides four additional (correlated) degrees of free- 
dom for each vector and axial- vector fit. 

5.1. Theoretical Uncertainties 

The estimate of theoretical uncertainties which 
is included in the fit takes into account uncer- 
tainties in the physical constants used and tiny 
contributions from quark mass corrections. In ad- 
dition, the effect of unknown higher order pertur- 
bative contributions is estimated by varying the 
leading unknown coefficient between zero and 



K4 ps 2\K 3 (K 3 /K2)\ ~ 50 (to be compared, e.g., 
with the experimental estimate K4 = 27 ± 5 ifyf ). 
The effect of the unknown fa coefficient of the j3- 
function is estimated in a similar way but turns 
out to be small. The uncertainty from the am- 
biguity of the renormalization scale \x is rather 
small; nevertheless it is considered here. We ad- 
ditionally include rough estimates of higher order 
non-perturbative operators which are not mea- 
sured here that lead to small contributions to the 
theoretical error. 

The dominant theoretical uncertainty for the 
prediction of R Tj v/A comes from the unknown 
coefficient. 

Results 

Computing the sum of the branching fractions 
from Table |] of the exclusive contributions to the 
vector and axial-vector spectral function yields 
the semileptonic widths 

R t V = 1.782 ± 0.018 (±0.016) , (19) 
R T , A = 1.711 ± 0.019 (±0.016) , (20) 
Rry+A+s = 3.649 ± 0.013 (±0.036) , (21) 

where the first errors are from experimental ori- 
gin while the second ones (in parentheses) rep- 
resent uncertainties in t he t heoretical predictions 

The value of © is 



as described in Section 5.1 



obtained via universality from the r branching ra- 
tios into e~v e v r and \x~VuV T and the r lifetime 
which are all taken from pT|,p| . The relative rise 
of the theoretical error in (|2l|) compared to jl9| ) 
and ( [2b"| ) comes from the additional uncertain- 
ties in the strange sector. It can be reduced to 
A th R T ^v+A+s — 0.034 when fitting simultane- 
ously the mass of the strange quark via R Ti s Jl5] |. 

Table || lists the measured values of the spec- 
tral moments ( |l8|) together with the experimental 
errors and the uncertainties from the theoretical 
predictions for the vector, axial- vector and vector 
plus axial- vector spectral functions. Note that 
all the moments obtained from the same spec- 
tral function are strongly correlated. The non- 
perturbative terms Jl7| ) are adjusted in the fit 
so that they do not contribute to the theoreti- 
cal errors. In particular, theoretical uncertainties 
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Table 3 

Fit results of a s (M T ) and the OPE non-perturbative contributions ( [i"7| ) from various 
input parameters: the V (A) fit uses R T y (Rt,a) and the corresponding moments. 
The (V + A) moments are fitted together with R t ,v+a+s- The combined fit exploits 
the information of R T y, R T ,A, Rt,v+a+s an d R T ,s as well as the vector and axial- 
vector moments. 



ALEPH 


Vector (V) 


Axial- Vector (A) 


V + A 


Combined Fit 
Vector Axial- Vector 


a s (M T ) 


0.354 ±0.020 


0.358 ±0.023 


0.348 ±0.017 


0.351 ±0.016 


S W 


(0.1±0.4)xl0 -3 


(-l.l±0.7)xl0 -3 


-(0.8±1.2)xl0 -3 


(-0.7±0.6)xl0 -3 


5 (6) 


0.029 ±0.004 


-0.028 ±0.004 


0.002 ±0.005 


0.023 ±0.003 -0.024 ±0.003 


5 (8) 


(-9.0±l.l)xl0 -3 


(7.7±1.0)xl0 -3 


(-1.0±1.0)xl0 -3 


(-8.8±1.0)xl0 -3 (8.3 ± 1.0) xlO -3 


xVd.o.f. 


0.4/1 


0.4/1 


0.1/1 


8.0/5 



Table 2 

Spectral Moments of vector (V), axial- vector 
(A) and vector plus axial-vector (V + A) r decay 
modes. The first error lines give the total exper- 
imental error including statistical and systematic 
uncertainties while the second error lines account 
for uncertainties in the theoretical prediction. 



ALEPH 


I = 


I = 1 


I = 2 


I = 3 


nU 

u v 


0.7159 


0.1689 


0.0532 


0.0227 


A c*P D U 


0.0034 
0.0037 


0.0006 
0.0035 


0.0007 
0.0004 


0.0006 
0.0002 


Da 


0.7205 


0.1471 


0.0639 


0.0303 


A ex P£) l( 

A th D\ l 


0.0033 
0.0037 


0.0009 
0.0031 


0.0005 
0.0005 


0.0004 
0.0003 


U V+A 


0.7177 


0.1581 


0.0585 


0.0265 




0.0022 
0.0037 


0.0006 
0.0033 


0.0004 
0.0004 


0.0004 
0.0002 



dominate the total error for the lower moments 
(I = 0, 1) due to the important error contribution 
from the missing K4 coefficient of the perturba- 
tive expansion. 

The results of the various fits of a s (M T ) and 
the non-perturbative terms of dimension D = 
4, 6, 8 are given in Table [|. The combined fit 
uses as input parameters R T .y, R T ,A, Rt.v+a+s, 
R T ,s and the vector and axial-vector moments. 
As mentioned above, the additional parameter 
R TjS = 0.1560 ±0.0064 || is used to adjust the 



strange quark mass in order to reduce the the- 
oretical uncertainty of the R T y+A+s prediction. 
The gain from the separation of vector and axial- 
vector channels compared to the inclusive V ± A 
fit becomes obvious in the adjustment of the lead- 
ing non-perturbative contributions of dimension 
D = 6 and D = 8, which cancel in the inclusive 
sum. The information for their accurate deter- 
mination comes mainly from the high I = 3, 4 
moments. Fig. || presents the total D = 6 contri- 
bution 5^ = C 6 (C(6))/M 6 to the vector and 
axial-vector hadronic widths as obtained from 
the combined fit. It is shown in comparison to 
recent experimental and theoretical constraints 
taken from |g|8||Q. 

Fig. H depicts the constraints of R T and the 
spectral moments on a s (M T ) when using stan- 
dard values for the non-perturbative contribu- 
tions H with the addition of large systematic un- 
certainties. They are in excellent agreement with 
the best a s (M T ) measurement obtained from the 
combined fit. The fact that the normalization of 
the respective spectra (R T ) and their pure shapes 
(spectral moments) independently give consistent 
results is a remarkable confirmation of the relia- 
bility of the a s measurement using r decays and 
the underlying SVZ approach (]l7|). 

The evolution @Q of the best a s (M T ), taken 
from the combined fit, to the Z boson mass yields 
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Figure 8. Non-perturbative contribution <5q from 
Eq. (0) to the r hadronic width of vector and 
axial-vector final states. The ellipse depicts the 
ALEPH result obtained in the combined fit. 
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Figure 9. Results for a s (M T ) using R T only, 
the moments only and the combined information 
from vector and axial-vector r decays. Addition- 
ally, the result from R T y + A+s only is plotted. 
The shaded band depicts a s (M T ) from the com- 
bined vector and axial-vector fit as given explic- 
itly in Table § 



the final result 

a s (M z ) = 0.1219 ± 0.0006 ± 0.0015 ± 0.0010 . 

The first error accounts for the experimental un- 
certainty, the second one gives the uncertainty of 
the theoretical prediction of R T and the spectral 
moments, while the last error stands for possi- 
ble ambiguities in the evolution due to uncertain- 
ties of the matching scales of the quark thresholds 
(0.0010) and effects associated with the trunca- 
tion of the RGE at next-to-next-to-leading order 
(0.0003). 

6. CONCLUSIONS 

In this talk, we presented new results on r 
hadronic spectral functions, measured by the 
ALEPH Collaboration at LEP, and their im- 
plications for QCD. With small inputs from 
unmeasured r decay modes, the vector and 
axial-vector contributions could entirely be sep- 
arated providing accurate experimental access 
to chiral sum rules and non-perturbative phe- 



nomcnological QCD. We used moments of the 
Laplace-transformed Das-Mathur-Okubo (DMO) 
sum rule for a precise determination of the DMO 
integral and, through it, of the electric polaris- 
ability of the pion. Finally, we used both vec- 
tor and axial-vector spectral functions for an ex- 
haustive study of a s and leading non-perturbative 
terms at the r mass scale. For the first time, 
the latter could be fixed with surprising preci- 
sion from r decays. A combined fit of all output 
variables yielded a s (M z ) = 0.1219 ± 0.0019 after 
the evolution from the r to the Z mass. 
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